{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "ACE Module\n",
    "==========\n",
    "\n",
    "This notebook demonstrates basic usage of the ``pyne.ace`` module."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import matplotlib.pyplot as plt\n",
    "import pyne.ace\n",
    "import requests\n",
    "url = \"https://www-nds.iaea.org/wolfram/w180/beta3/W180.ace\"\n",
    "r = requests.get(url)\n",
    "with open(\"W180.ace\", \"wb\") as outfile:\n",
    "    outfile.write(r.content)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The main class in ``pyne.ace`` is called ``Library``. It is instantiated using the name of an ACE file, in this case one distributed with MCNP."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "lib = pyne.ace.Library('W180.ace')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One can choose to read all tables in the file or selectively read a subset by specifying an argument to the ``read`` method."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "lib.read('74180.21c')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "After the call to ``read()``, the Library instance will have a dictionary called ``tables``."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'74180.21c': <ACE Continuous-E Neutron Table: 74180.21c>}"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lib.tables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "w180 = lib.tables['74180.21c']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Once a table is selected, we can inspect, e.g., the energy grid and the total cross section."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1.00000e-11, 1.03125e-11, 1.06250e-11, ..., 5.00000e+01,\n",
       "       5.50000e+01, 6.00000e+01])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "w180.energy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1544.50481   , 1520.95229   , 1498.44763   , ...,    4.17220824,\n",
       "          4.26476628,    4.3440783 ])"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "w180.sigma_t"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To get data on a reaction, such as fission or $(n,2n)$, there is an attribute called ``reactions``."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "OrderedDict([(2, <ACE Reaction: MT=2 (z,elastic)>),\n",
       "             (5, <ACE Reaction: MT=5 >),\n",
       "             (16, <ACE Reaction: MT=16 >),\n",
       "             (17, <ACE Reaction: MT=17 >),\n",
       "             (22, <ACE Reaction: MT=22 >),\n",
       "             (24, <ACE Reaction: MT=24 >),\n",
       "             (28, <ACE Reaction: MT=28 >),\n",
       "             (41, <ACE Reaction: MT=41 >),\n",
       "             (51, <ACE Reaction: MT=51 >),\n",
       "             (52, <ACE Reaction: MT=52 >),\n",
       "             (53, <ACE Reaction: MT=53 >),\n",
       "             (54, <ACE Reaction: MT=54 >),\n",
       "             (55, <ACE Reaction: MT=55 >),\n",
       "             (56, <ACE Reaction: MT=56 >),\n",
       "             (57, <ACE Reaction: MT=57 >),\n",
       "             (58, <ACE Reaction: MT=58 >),\n",
       "             (59, <ACE Reaction: MT=59 >),\n",
       "             (60, <ACE Reaction: MT=60 >),\n",
       "             (61, <ACE Reaction: MT=61 >),\n",
       "             (62, <ACE Reaction: MT=62 >),\n",
       "             (63, <ACE Reaction: MT=63 >),\n",
       "             (91, <ACE Reaction: MT=91 >),\n",
       "             (102, <ACE Reaction: MT=102 >),\n",
       "             (103, <ACE Reaction: MT=103 >),\n",
       "             (107, <ACE Reaction: MT=107 >),\n",
       "             (203, <ACE Reaction: MT=203 >),\n",
       "             (207, <ACE Reaction: MT=207 >),\n",
       "             (444, <ACE Reaction: MT=444 >),\n",
       "             (600, <ACE Reaction: MT=600 >),\n",
       "             (601, <ACE Reaction: MT=601 >),\n",
       "             (602, <ACE Reaction: MT=602 >),\n",
       "             (603, <ACE Reaction: MT=603 >),\n",
       "             (604, <ACE Reaction: MT=604 >),\n",
       "             (605, <ACE Reaction: MT=605 >),\n",
       "             (606, <ACE Reaction: MT=606 >),\n",
       "             (607, <ACE Reaction: MT=607 >),\n",
       "             (608, <ACE Reaction: MT=608 >),\n",
       "             (609, <ACE Reaction: MT=609 >),\n",
       "             (610, <ACE Reaction: MT=610 >),\n",
       "             (611, <ACE Reaction: MT=611 >),\n",
       "             (612, <ACE Reaction: MT=612 >),\n",
       "             (613, <ACE Reaction: MT=613 >),\n",
       "             (614, <ACE Reaction: MT=614 >),\n",
       "             (615, <ACE Reaction: MT=615 >),\n",
       "             (616, <ACE Reaction: MT=616 >),\n",
       "             (617, <ACE Reaction: MT=617 >),\n",
       "             (618, <ACE Reaction: MT=618 >),\n",
       "             (619, <ACE Reaction: MT=619 >),\n",
       "             (620, <ACE Reaction: MT=620 >),\n",
       "             (621, <ACE Reaction: MT=621 >),\n",
       "             (622, <ACE Reaction: MT=622 >),\n",
       "             (623, <ACE Reaction: MT=623 >),\n",
       "             (624, <ACE Reaction: MT=624 >),\n",
       "             (649, <ACE Reaction: MT=649 >),\n",
       "             (800, <ACE Reaction: MT=800 >),\n",
       "             (801, <ACE Reaction: MT=801 >),\n",
       "             (802, <ACE Reaction: MT=802 >),\n",
       "             (803, <ACE Reaction: MT=803 >),\n",
       "             (804, <ACE Reaction: MT=804 >),\n",
       "             (805, <ACE Reaction: MT=805 >),\n",
       "             (806, <ACE Reaction: MT=806 >),\n",
       "             (807, <ACE Reaction: MT=807 >),\n",
       "             (808, <ACE Reaction: MT=808 >),\n",
       "             (809, <ACE Reaction: MT=809 >),\n",
       "             (810, <ACE Reaction: MT=810 >),\n",
       "             (811, <ACE Reaction: MT=811 >),\n",
       "             (812, <ACE Reaction: MT=812 >),\n",
       "             (813, <ACE Reaction: MT=813 >),\n",
       "             (814, <ACE Reaction: MT=814 >),\n",
       "             (815, <ACE Reaction: MT=815 >),\n",
       "             (816, <ACE Reaction: MT=816 >),\n",
       "             (817, <ACE Reaction: MT=817 >),\n",
       "             (818, <ACE Reaction: MT=818 >),\n",
       "             (819, <ACE Reaction: MT=819 >),\n",
       "             (820, <ACE Reaction: MT=820 >),\n",
       "             (821, <ACE Reaction: MT=821 >),\n",
       "             (822, <ACE Reaction: MT=822 >),\n",
       "             (823, <ACE Reaction: MT=823 >),\n",
       "             (824, <ACE Reaction: MT=824 >),\n",
       "             (825, <ACE Reaction: MT=825 >),\n",
       "             (826, <ACE Reaction: MT=826 >),\n",
       "             (827, <ACE Reaction: MT=827 >),\n",
       "             (828, <ACE Reaction: MT=828 >),\n",
       "             (829, <ACE Reaction: MT=829 >),\n",
       "             (830, <ACE Reaction: MT=830 >),\n",
       "             (831, <ACE Reaction: MT=831 >),\n",
       "             (832, <ACE Reaction: MT=832 >),\n",
       "             (833, <ACE Reaction: MT=833 >),\n",
       "             (834, <ACE Reaction: MT=834 >),\n",
       "             (849, <ACE Reaction: MT=849 >)])"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "w180.reactions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "elastic = w180.reactions[2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "An instance of a Reaction contains the reaction cross section and any angular or energy distributions that may be present."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([43.83981, 43.20029, 42.58963, ...,  1.9508 ,  2.0991 ,  2.2343 ])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "elastic.sigma"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With the energy grid (stored on the table), and the cross section (stored on the reaction), one can generate plots of the cross section."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7f5e13830ef0>]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAlE0lEQVR4nO3deXyU1f328c83k4SQkEQIYU9YDAQRBDUgiMhSRCwi7krVVosirdhitVYf+3vap9a61KVqXYqKqFWo9acWLEotoggii+ICsogsGpBdWYVs5/kjmRBiEiazZO7MXO/XC2Hu3HPPd2JyzZlzzpxjzjlERCT2JUS7ABERaRgKfBGROKHAFxGJEwp8EZE4ocAXEYkTCnwRkTiRGO0C6tKyZUvXqVOnaJchItJofPDBBzucc9k1fc2TgW9mo4HReXl5LF26NNrliIg0Gma2sbavebJLxzk30zk3PjMzM9qliIjEDE8GvoiIhJ8CX0QkTijwRUTihAJfRCROKPBFROKEAl/kKNZs3UtJaVm0yxAJmScD38xGm9nk3bt3R7sUiXNrt+1jxAPzuO/NNdEuRSRkngx8zcMXr9i25yAAy778JsqViITOk4Ev4hX+/eAMi2odIuGgwBcJgCnvJQYo8EUCoK2fJRYo8EXq4G/YO5T40vgp8EXqUpH4obTw//7+RrbsPhieekRCoMAXqYN/sDbYvN+25yC/fXU5V01dEr6iRIKkwBepgx3u0wlKSVn5Hb89UBSegkRC0GCBb2bHmdnjZvaSmf2soR5XJBThmpyjQV/xgpAC38ymmNk2M1te7fhIM1ttZmvN7BYA59xK59wE4GJgYCiPK9LQgh201XRO8ZJQW/hTgZFVD5iZD3gEOAvoAYw1sx4VXzsH+DcwK8THFWkQVpHYwbbQ1bIXLwkp8J1z84Bd1Q73A9Y659Y554qA6cCYivNnOOfOAi4L5XFFGopa6BJLIrGJeXvgqyq3C4FTzGwIcD7QhDpa+GY2HhgPkJubG4HyROov2Ia6XjDESyIR+DVyzr0NvB3AeZOByQAFBQV6QyxRVTlJR30zEgMiMUtnE5BT5XaHimMB0/LI4hVqoUssiUTgLwG6mllnM0sGLgVm1OcCWh5ZvEbte4kFoU7LnAYsBPLNrNDMxjnnSoCJwGxgJfCic25FPa+rFr54hJr4EjtC6sN3zo2t5fgsQph66ZybCcwsKCi4JthriIRTqF34WnxNvEBLK4jUwd+HH/QsHb1DEA/xZOCrS0e8QnEtscSTga9BW/EcTcuUGODJwBfxGsW9xAJPBr66dMQrwhX0eoMgXuDJwFeXjniFglpiiScDX8RrQg1+fWJXvMCTga8uHfEOV/Hf0BJf7xTECzwZ+OrSEa9Qy15iiScDX8RrtAGKxAIFvkgd/Hmt4JZYoMAXqYO6dCSWeDLwNWgrXuHf+CTUBr7eIIgXeDLwNWgrXhPsjldq4IuXeDLwRbxCLXuJJQp8kTposFZiiQJfJADBBr+6dMRLPBn4GrQVr3Bh+qStiBd4MvA1aCueoZyXGOLJwBfxCuW9xBIFvkgANHgrsUCBL1IHf9CHPD1TLxjiAQp8kTposFZiiQJfpA6VLXw10SUGKPBFAhB03GsivniIJwNf8/DFK9Sul1jiycDXPHzxCheuUVsRD/Bk4It4jfJeYoECX6QOCnqJJQp8kbqEbZaOXjok+hT4InXQPHyJJQp8kQAEHft6vRAPUeCL1OHwB6/Ccx2RaFLgi9ShvkF9oKiE977Y8f3rhKkekVAo8EUCEGhf/k3//JgfPbGITd9+d+T91cQXD0hsqAcys3OBUUAG8JRz7j8N9dgiwfLHdKB5vXrLXgAOHCo54v4iXhBSC9/MppjZNjNbXu34SDNbbWZrzewWAOfcq865a4AJwCWhPK5IQ6lvy9ys5sVzFPziBaF26UwFRlY9YGY+4BHgLKAHMNbMelQ55bcVXxfxvPq28GtTVqbIl+gLKfCdc/OAXdUO9wPWOufWOeeKgOnAGCt3N/C6c+7D2q5pZuPNbKmZLd2+fXso5YlETfV4V9yLF0Ri0LY98FWV24UVx64HhgMXmtmE2u7snJvsnCtwzhVkZ2dHoDyRwNW3ZV+9Q0djteIlDTZo65x7CHgokHPNbDQwOi8vL7JFiRxVeWLXty//e6cr+MUDItHC3wTkVLndoeJYwLQ8snhNqHm991AJL31QGJZaRIIVicBfAnQ1s85mlgxcCsyIwOOIRFw4u2Ru+ufHlGrwVqIo1GmZ04CFQL6ZFZrZOOdcCTARmA2sBF50zq2o53W145V4Qn1n6fhnZdb2QS19AEuiKaQ+fOfc2FqOzwJmhXDdmcDMgoKCa4K9hkg4hLyGjjrvxUM8ubSCWvjiNQpuiQWeDHwN2opXuMpZOoGdbxUTM2s7v6TMqVtHoqbBpmWKNEb1nodf88oKlfr84T+0zWzKyJ5t+GHPtvRsn1Hrcgwi4ebJFr66dMQrQt7YsNoFDhaXUVRSxuR56xj91/mMeWQBrywrpLi0LMRHEjk6Twa+unTEa+ob/HW9M5g4LI+ltw3n9jHHs/9QCTf842POuP8d3li+JaQaRY7Gk4Ev4hX+/vZw97s3T0vmigGdePOGwTzx4wKSExOY8PcPuO75D9m571BYH0vET4EvEoCSCH1gKiHBOKNHa2b9YhC/PjOf/3y2hVEPzWfFZnVnSvh5MvDVhy9eU1xSvz72ytk9AZ6f6EvguqF5vPLzgZjBRY8v5O3V2+pZpUjdPBn4ofbh3zt7Nf/6qF7L94jUyN+TU1waWHSHOuOmZ/tMXr1uIJ2y0rj6maW8ukw/xxI+ngz8UBSVlLFkwy5+Of0jHpm7VnOeJST+lnpxWcPNommdkcL0a/tT0Kk5k/7xEVc/s5TPNu9psMeX2BVzgZ+cmMCz4/oxpk87/jx7Nbe+/KmmvEnQ/O2F+rYbQm1nZKQkMfWqftx4RjcWrd/JqIff5fbXPuO7otLQLixxzZOBH2offpNEH3+5pA8Th+YxfclXjHtmKfsqNpUWiaRwfoQqJcnH9T/oyvybh3HZKbk8NX89Ix+cxwcbvwnjo0g88WTgh2Mevplx05n53HV+Lxas3cFFjy9ky+6DYaxS4kHIi6eFoUsxMzWJP57bi+nj+1NS6rjo8ff406yV7D1YHPK1Jb54MvDD6dJ+uUy5si9f7tzPeY8uYNUW9YVK4Oob15FcJaF/lyxenzSIi07OYfK8dQy77x1mfrxZ41QSsJgPfIDB3bL554RTcQ4ufGwh736uzdGlccpISeLuC0/g1esG0jYzheunLeOaZ5fy9e7vol2aNAJxEfgAPdpl8Mp1p9KheVOuenoJLyz6MtolSSNQtfVcn5Z0pBvdfXKO4ZWfD+S3o45j/todnHH/PJ5buEE7akmdPBn4kfrgVdvMpvxzwgAG5rXk/7zyKb/713LN4JE6VY3PogB+VhIq+nTKKpdkiERV5XwJxtWDuvCfSYPpk3MM//OvFYx+eD6L1u2M3INKo+bJwI/k4mnpKUlMubIv1wzqzDMLN/KTKYv59kBR2B9HYkSVwC4K4NO2CZVbHDac3KxUnhvXj4fHnsi3B4q4ZPL7XPfChxR+c6ABq5DGwJOBH2m+BOO2UT2496LeLN3wDWMeWcDnW/dGuyzxuO+Kjz4H3qq18BuKmTG6dzvm3DiEScO7MmflVn5w3zvc/+Yazd2XSnEZ+H4XntyBaeP7s/9QKec9+h5vrdoa7ZLEY6pubbjv4NE/y1HZwo/SzJmmyT4mDe/GnBuHMOL4Njw053OG3fc2MzSbR4jzwAc4uWNzZkwcSMesVMY9s5QH3lyjgS+pVDUjA/nw3uE+/EhVFJj2xzTl4bEn8uK1A2iRlswvpi3joscX8mmhFiSMZ3Ef+ADtjmnKSxNO5bw+7Xlwzudc+fRirUku3xNYC78i8KOd+BX6dW7BjImncdf5vdiwcz/nPDKf37z0Cdv36uc7HinwKzRN9nHfxb258/xeLFq/i1EPzWfphl3RLkuirGps7w2ghe//4JVH8h4oH7O6tF8ub900hKtP68z/fljI0HvfZvK8LwIaiJbYocCvwswY2y+Xl392Kk2SErhk8vv87Z0vPNNak4Z3RJdOAC18i3Iffl0yUpK4bVQPZt9wOv06t+BPs1Zx5l/mMWflVk/WK+HnycCP9gYoPdtnMvP60xjRozV3vr6KHz35vqa4xamqg7Z7Ali7xit9+HU5NrsZU67sy9NX9cUMxj2zlMufWsTyTerfj3WeDHwvbGKekZLEo5edxD0XnsCnhbs56y/v8vKHhWoJxbFA+r0b8oNXoRqa34rZk07nd6N78NnmPYz+63x+9eJHbP5WyzTEKk8GvleYGRcX5PDGpNPp3jadX734Mdc+94HWLYkj/sBOTDC27gkg8BOiMw8/WEm+BK4a2Jm3fz2U8ad34bVPvmbovW/zx9c+Y9terS4baxT4Achpkcr08QO49azuzPt8O8Pve4en5q+nRMsyxDx/bLfOSAkoAP2LZTaSvK+U2TSJW886jrduHMyoE9oyZcF6Bt09l9sV/DFFgR8gX4Jx7eBjefOGwfTt3ILbX/uMcx9doHVL4kSbzBS+DmA/hcNLKzSyxK/QoXkq91/chzk3DuHsE9ox9b0NDLp7Ln+Y+Rnb9ij4GzsFfj3ltEjl6Sv78siPTmLH3vJ1S8ZNXcIaLc0Qmyqa6p1bprFx5/6jLrbn78P3n9ZYg79zyzTuu7g3c341mNG92/HMwg0Mumcu/2/mCgV/I6bAD4KZMeqEtsy9aQg3j8xn8fpdjPzLPCZNX8bKr7XBSizxx3X3NukUlzrW79hf5/lV19JZ9uU3jX6XtU4t07j3ot68deNgzundjmcXblTwN2KJ0S6gMWua7OPnQ/IY2zeXx975guff38irH21maH42Pz2tMwOPbVk5iCeNk78v/ri2GQCs/HoP3Vqn13r+4bV04LxH34t0eQ2mY1Yaf76oNxOH5fHI3LU8u3Ajzy/6kvNPbM/VgzqT16r274l4h1r4YdA8LZn/88PjeO+WH3DTiG58UribK55azOl/nstDcz7nq12aw9/YdW3djIyUROZ/vqPO8/xdOnVN3z0UwKqbXtUxK417LuzN3BuHcNHJHXhl2SaG3z+Pq55ezHtrd2jasscp8MMoMzWJicO6suCWYTw89kQ6ZaVx/5trGHTPXM568F0eeHMNyzft1uJsjYg/wBITEhjULZu312yv8/9fUmL5r9ShOpYsuGPWyvAWGQW5WanccV4v3rtlGDcML2/k/OjJRfzwoflMX/yllmT2qAbr0jGzLsBtQKZz7sKGetxoSEnyMbp3O0b3bsdXuw4we8UWZq/YwkNvfc6Dcz4nPSWRvp1a0K9zC3q1z6R7m3SymjWJdtlSA3+0G3B2r7b8+5OvmbNyKyOOb1Pj+c2alP9K1bWyZnFp7LzgZzVrwi+Hd+XawV14ddkmpr63gVte/pS73ljFJQU5XN6/IzktUqNdplQIKfDNbApwNrDNOdezyvGRwIOAD3jSOXeXc24dMM7MXgrlMRubnBapXD2oC1cP6sL2vYeYv3Y7i9fvYtH6Xby1alvlea3Sm5DfJp2cFqnkNE+lQ/OmdGjelNYZKbRISyYlyRfFZxG//D0UZnBGj9a0P6YpD/z3c4Z1b0Wi7/tvkNNTjh74sSglycel/XK5pG8Oi9bv4tmFG3hy/nomv7uOYfmtuHxARwZ3zdaYVpSF2sKfCvwVeNZ/wMx8wCPAGUAhsMTMZjjnPgvxsRq97PQmnHdiB847sQMAO/cdYuXXe1m1ZQ8rv97Lmq17Wb7pa7458P01W1KTfbRIS6ZFWjKZTZNIS04ktYnvyL+TfaQk+UjyGUm+BBJ9CSRX+bf/eJIvgcQEI9FnJCYYvoSEir+tyt8J+Hx2xHH/DJR4ZBiJvgT+5+zjmPD3D3ngv2v49Zndv3deWnJF4Aew0FosMjP6d8mif5csNn/7HdMWf8m0xV8x5+kl5LZI5fL+uVx0cg7N05KjXWpcCinwnXPzzKxTtcP9gLUVLXrMbDowBoj7wK8uq1kTTuvahNO6tjzi+L5DJWz65ju+2nWA7fsOsWt/UeWfnfuL2H2giC27D3KgqJQDRSXsLyptkGVuE6y8L7vyhcFX7QWiyguGr+IFpeYXkyrHfTUfT/QZ6SlJZKQkktE0iYyUJDKaJpKRkkSLtGSy05uQVEMLO9yqd76ceXwbLu2bwyNzv6BJoo/rh+Ud8UKY6Cv/9879tS/DkNeqWSRK9Zx2xzTlxhH5XD+sK7NXbOG5hRv506xV3PefNYzu3Y4r+nekd84x0S4zrkSiD7898FWV24XAKWaWBdwBnGhmtzrn7qzpzmY2HhgPkJubG4HyvK9Zk0Ty26ST3ybwqW4lpWUcKC7lYFEpxWWOktIyikvLKCpxlJSV/7u41FFcWkZJqaOotIyyMkdJmaO08u+yw7dLazl+xNdrOO4/v7T68fIaviuufv2yGu7vKCopq7NbxAxapCbTKiOF1hlNaJvZlM4tU+mYlUbnlmnktkgNSzdY5awT8z+ucfu5PSkqKeP+N9fwSeFu/nR+T1qlpxxxv7Xb9tV4vaev7MvAvJY1fi1WJScmVI5prdqyh+cWbuSVZZt46YNCTuiQyRX9OzK6dzt1WzaABhu0dc7tBCYEcN5kYDJAQUFB7IxuRViiL4EMXwIZKUnRLiVsSssc+w6VsOe7YvYcLGbPdyXs/q6YXfuL2Lb3IFv3HGJ7xd+fFO5m1/6iyvuaQeesNHq2z6RX+0x6ts/khA6ZpDUJ7ke+am9Wki+B+y7uTY92GdwzezUjHpjHn87rxQ97ta08Z9mX39Z4naHdWwX1+LGie5sM7jivF7ec1Z1Xlm3iuYUb+fVLn3DHrJVcXJDDZafk0jErLdplxqxIBP4mIKfK7Q4VxwJmZqOB0Xl5eeGsSxoZX4KR2TSJzKaBvYjtPlDMhp372bBzP+t37OezzXtYumEXMz7eDJSveHlSx+ac3rUlQ/JbcXy7jKDHJcyMqwd1YUh+K2588SN+/vyH/HhAR7IrZluVaOptndJTkvjxgE5c0b8j76/bxd/f38iU+euZPG8dP+jeiknDu9GrQ/SWR49VFuoHJSr68F/zz9Ixs0RgDfADyoN+CfAj59yK+l67oKDALV26NKT6RHbsO8Snm3azeP0u3v18O8s3lS9/0SU7jQtO6sDlp3QkM7XmF5Un5q3jjlkr+fT3I0iv5d1TcWkZd7++iifnryc5MYGikjJaZzSpcTnlDXeNCt8TizFb9xzkhUVf8szCDXx7oJjhx7Vm0vCu9Gyv4K8PM/vAOVdQ09dCGvUys2nAQiDfzArNbJxzrgSYCMwGVgIv1jfso73jlcSWls2aMDS/Fb8Z2Z3Xrh/EB78dzl3n96JlWhP+PHs1p939Fg+8uYZDJd//sJB/8bO63gkk+RK4bdRxjOrVtnLw/KKTc2o9X2rWOiOFG87oxrs3D+XGM7qxeP1Ozn54Ptc+t1RrVIVJSIHvnBvrnGvrnEtyznVwzj1VcXyWc66bc+5Y59wdQVw36jteSezKataES/vl8uKEAbwxaRCndW3Jg3M+Z8xfF7ChlsXRjtbxY2ZcMaBj5e3Tu2WHseL4kp6SxPU/6Mq7vxnGpOFdeW/tTs568F2ue/5DrUobIk8uraAWvjSU7m0yeOzyk5lyZQFb9xzk4r8tPCL069Pj2atK10PbzJQ6zpRAZDZNYtLwbsz/zTCuH5bHO2u2c+Zf5nH9tGW1zoKSunky8NXCl4Y2rHtr/nHtAIpKy7juhQ8ru3cql1YIYGy36gwgTTEMn8zUJG4ckc+7Nw/lZ4OPLV/a4oF3uPHFj7UwYT15MvBFoqFb63TuvuAEVmzew/TFXx3xNTtqp86RUpL0qxVuzdOSuXlkd969eSjjTuvMa59sZui9b3PbK59qn+kAefKnUl06Ei0jerTmlM4teOztLygtc0HvTdsQnwKOV1nNmnDbqB688+uhjO2Xy4tLv2Lwn9/mDzM/Y/veo280H888+VOpLh2JFjPj8v4d2bLnIIvW76wyS6d+10nUImER1yYzhdvP7clbNw7h3D7l2zCefs9c7n5jFd8eKDr6BeKQJwNfJJqGH9eaJokJzFm5LegWfjwvNNfQclqkcs+FvXnzhtMZcXxrHn/nCwbdPZcH//s5ew9+fyHCeKbAF6mmabKP3jnHsHTjN5XHgsnvd28eSqt07XPQULpkN+PBS0/kjV+ezql5WTzw3/LNhx5/5wsFfwVPBr768CXaTu7YnBWbdoe0CmlOi1SyFfgNLr9NOn+7ooCZE0+jT84x3PX6Kgbc+Ra/n7HiqJvQxzpPBr768CXaurZqRkmZq5z2V99ZOv6zU5M1PTNaenXIZOpV/fjXdQM5o0drnl+0kaH3vs1Ppy6J2+D3ZOCLRFunluUrNq6rCIZgu+Q1Wyf6euccwwOX9GFBxf67SzbsYuRf5vHEvHVxt7+0fhpFatC5YonejTtDawn6A//mkflMvapvyHVJ8Fqlp/DL4V35768GM6hrS+6YtZIfT1kUV/37ngx89eFLtB2TmkRyYkLldpP1beD73xEkVeyA1aF5KkPy43stfK9onZHCEz8u4O4LerFo3S4ueOw9nnx3HZ8W7o75Fn+DbYBSH865mcDMgoKCa6Jdi8QnM6NlWjKbdx+svB0MfwtfkzS9xcy4pG8ubTKb8vsZK/jjv1cCkN4kkYJOzRmY15ILTuoQc3vvejLwRbwgq1mTw4Ef5DX8gV8W4r4TEhmDu2Uz96YhbNld/kG7Ret3sWjdTub+eyUPvLmGX5+ZzxUDOuGLkQ/SKfBFatGyWfCtO/87An/gN8Qm8xK8NpkpjOnTnjF92gOwZute/vjvlfx+5me8tXo7j112UtDbY3qJJ/vwRbwgo8rWisHP0im/Y3GpWviNSbfW6TxzVV/uOK8nC9bu4EdPLmLfoZJolxUyBb5ILVKTD7fogu3DT6joClCXTuNjZlx2Skceu+wklm/azS+nLWv0g7qeDHzN0hEvaNYk9A9N+bt+Q907WqJnxPFt+P05xzNn1TbuqBjcDcTu74qZu3obC9buqHH7zGjwZKeUZumIF1Rt4Qcrwfwt/JAvJVF0Rf+OrN++nykL1tMlO43L+3es9dyyMsdj73zBo3PXsr+oPOhzW6Ry9wUnMODYrIYquUaebOGLeEGzMAzS+TuC1MJv/G4bdRzDurfidzNW8N7aHTWec6CohJ89/wF/nr2a07q2ZNo1/fnbFSeT6DN+8vRiFtRyv4aiwBepRWoYunT8ff+K+8bPl2A8eGkfurRM42fPf8jqLUduqL752++4+G8LefOzrfzfs3vw+OUnM+DYLM48vg3/O+FUOmWl8otpy9ixL3qbtCjwRWqRFoYuHf9Yr7p0YkN6ShJP/aQvKUkJXDJ5IbM+/Zr9h0p4Y/nXnPPX+WzYcYAnf1LAT0/rfMRAf/O0ZB4eexJ7D5Zw2yuf1vmOb8/B4oi9E1Dgi9SiSWLovx7+Pnx16cSO3KxU/nntqbTNbMrPn/+Q4383mwl//5AWacm8et2pDOveusb75bdJ56YzuzF7xVZe/WhTjec457jlfz/hyqcXR2SfXk8O2op4QThWukyobOEr8GNJblYqMyYOZM7KrXyxfT/HZqcx/LjWJB7lZ2bcaV14ffkW/jDzM07vmk1WsyP3S/j7+xuZ9ekWbjmrO20zm4a9brXwRWqRFIYWvmmWTsxK8iUwsmdbrhuax8iebY8a9lA+DnDX+Sew71BJ5fo9fis27+b211YyJD+b8YO6RKRmTwa+5uGLF/g/JRsKq5yHH/KlJEbkt0lnwuBjeWXZJt5atRWAfYdKmPjCMpqnJXH/xX0qP7AXbp4MfO14JV6QHESXzu9H9yCrygqLh+fhK/HlsOuG5tG9TTrXPb+Mm/75Mec/uoCNO/fz0KUn0iKCK3R6MvBFvCCYPvwrB3bmg/85o/K2PmkrNUlJ8vHMT/txVq82zF6xhW8OFPP0Vf04pUtkP5ilQVuRWoRj0Na/F6768KW61hkp3H9xH4pLy/CZRawbpyoFvkgtkhPVhy+R15D7HqtLR6QWyb7wfdJWffjiBQp8kVokhaGFX9mHH/KVREKnwBepRTj78DVoK16gwBepRWIYBtH0SVvxEgW+SC3CMWsi2K0RRSKhwWbpmFka8ChQBLztnHu+oR5bJBgJYUxrNfDFC0Jq4ZvZFDPbZmbLqx0faWarzWytmd1Scfh84CXn3DXAOaE8rkhD8IUh8IPdC1ckEkLt0pkKjKx6wMx8wCPAWUAPYKyZ9QA6AF9VnOaNDR5F6hDOrFYDX7wgpMB3zs0DdlU73A9Y65xb55wrAqYDY4BCykO/zsc1s/FmttTMlm7fvj2U8kRC4muATz6KNKRIDNq253BLHsqDvj3wMnCBmT0GzKztzs65yc65AudcQXZ2dgTKEwlMOPvwRbygwQZtnXP7gasCOdfMRgOj8/LyIluUSB3UwJdYE4kW/iYgp8rtDhXHAqblkcULNOAqsSYSgb8E6Gpmnc0sGbgUmFGfC2gDFIk1mpYpXhDqtMxpwEIg38wKzWycc64EmAjMBlYCLzrnVtTnumrhS6zQmwTxkpD68J1zY2s5PguYFcq1RUQkvDy5tIK6dCTWOM3EFw/wZOCrS0dihX+1TBEv8GTgi8QcNfDFAzwZ+OrSEREJP08Gvrp0RETCz5OBLxIrNC1TvMSTga8uHYk16sIXL/Bk4KtLR2KFGvjiJZ4MfBERCT8FvkgDcFpMRzxAgS8SQRq0FS/xZOBr0FZEJPw8GfgatBURCT9PBr5IrFEXvniBAl8kgrR4mniJAl9EJE54MvA1aCuxRj064gWeDHwN2kqs0LRM8RJPBr5IrNGgrXiBAl9EJE4o8EVE4oQCX0QkTijwRRqA0zwd8QAFvkgEmabpiId4MvA1D19EJPw8Gfiahy+xRtMyxQs8GfgisUIdOuIlCnyRCBpwbBYAI45vHeVKRCAx2gWIxLLj2maw4a5R0S5DBFALX0QkbijwRUTihAJfRCROKPBFROKEAl9EJE40WOCbWRcze8rMXmqoxxQRkcMCCnwzm2Jm28xsebXjI81stZmtNbNb6rqGc26dc25cKMWKiEjwAp2HPxX4K/Cs/4CZ+YBHgDOAQmCJmc0AfMCd1e7/U+fctpCrFRGRoAUU+M65eWbWqdrhfsBa59w6ADObDoxxzt0JnB1sQWY2HhgPkJubG+xlRESkmlD68NsDX1W5XVhxrEZmlmVmjwMnmtmttZ3nnJvsnCtwzhVkZ2eHUJ6IiFTVYEsrOOd2AhMCOdfMRgOj8/LyIluUiEgcCaWFvwnIqXK7Q8WxkGl5ZBGR8Asl8JcAXc2ss5klA5cCM8JTloiIhFug0zKnAQuBfDMrNLNxzrkSYCIwG1gJvOicWxGOorTjlYhI+AU6S2dsLcdnAbPCWlH5dWcCMwsKCq4J97VFROKVJ5dWUAtfRCT8PBn4GrQVr/jFsDyuPLVTtMsQCQvteCVSh1+NyI92CSJh48kWvrp0RETCz5OBry4dEZHw82Tgi4hI+CnwRUTihCcDX334IiLh58nAVx++iEj4eTLwRUQk/BT4IiJxwpxz0a6hVma2HdgY7Tpq0RLYEe0igtSYawfVH02NuXaIj/o7Oudq3D3K04HvZWa21DlXEO06gtGYawfVH02NuXZQ/erSERGJEwp8EZE4ocAP3uRoFxCCxlw7qP5oasy1Q5zXrz58EZE4oRa+iEicUOCLiMQJBb6ISJxQ4IeJmXUxs6fM7KUqx9LM7Bkze8LMLotmfYEys1wze9XMppjZLdGup77MLMHM7jCzh83sJ9Gup74qfmaWmtnZ0a6lvszs3Iqf9X+Y2Yho1xOIxvg76hfM91uBD1SE2zYzW17t+EgzW21ma48Wfs65dc65cdUOnw+85Jy7BjgnzGV/TzieB9CL8pp/CpwYsWJrEKb6xwAdgGKgMFK1Vhem2gF+A7wYmSprF6bfgVcrftYnAJdEst661PO5NOjv6NHUp/agvt/Oubj/A5wOnAQsr3LMB3wBdAGSgY+BHpQH4mvV/rSqcr+Xqvz7VqBPxb9faAzPA8gC5gJvAVc1tv8PwC3AtdX/XzSS2s8ALgWuBM5ubN/7Kve7DzipIesP4bk06O9oOGsP5vutTcwB59w8M+tU7XA/YK1zbh2AmU0Hxjjn7gQCfbtdSHlr8yMa4N1UOJ6Hmd0E/K7iWi8BT0e47Ephqr8QKKq4WRrBco8QptqHAGmUB9F3ZjbLOVcWybr9wlS/AXcBrzvnPoxwybWqz3OhgX9Hj6Y+tZvZSur5/Y76E/Sw9sBXVW4XVhyrkZllmdnjwIlmdmvF4ZeBC8zsMWBmxCqtW72eB/AG8IuK57IhgnUFqr71vwycaWYPA/MiWVgA6lW7c+4259wk4AXgiYYK+zrU93t/PTAcuNDMJkSysCDU9ly88Dt6NLXVXu/vt1r4YeKc20l5X1rVY/uBq6JTUXCcc8uBC6NdR7CccweA6mMpjYpzbmq0awiGc+4h4KFo11EfjfF31C+Y77da+LXbBORUud2h4lhj09ifR2OuvzHXDo2//qoa83MJW+0K/NotAbqaWWczS6Z8MG1GlGsKRmN/Ho25/sZcOzT++qtqzM8lfLVHe1TaC3+AacDXHJ7KN67i+A+BNZSPkN8W7Tpj/Xk05vobc+2xUH+sPJdI167F00RE4oS6dERE4oQCX0QkTijwRUTihAJfRCROKPBFROKEAl9EJE4o8EVE4oQCX0QkTijwRUTixP8HBZ6U0aYL3cYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.loglog(w180.energy, elastic.sigma)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
